GeV Photons from Ultra High Energy Cosmic Rays accelerated in Gamma Ray Bursts 
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7-ray bursts are produced by the dissipation of the kinetic energy of a highly relativistic fireball, 
via the formation of a collisionless shock. When this happens, Ultra High Energy Cosmic Rays 
up to ~ 10 20 eV are produced. I show in this paper that these particles produce, via synchrotron 
emission as they cross the acceleration region, photons up to 300 GeV which carry away a small, 
~ 0.01, but non-negligible fraction of the total burst energy. I show that, when the shock occurs 
with the interstellar medium, the optical depth to photon-photon scattering, which might cause 
energy degradation of the photons, is small. The burst thusly produced would be detected at Earth 
simultaneoulsy with the parent 7-ray burst, although its duration may differ significantly from that 
of the lower energy photons. The expected fluences, ~ 10~ 5 — 10 -6 erg crn~ 2 , are well within the 
range of planned detectors. A new explanation for the exceptional burst GB 940217 is discussed. 



PACS numbers: 98.70.Rz, 98.70.Sa 

It can be argued compellingly that a model for ex- 
tragalactic 7-ray bursts (GRBs) must involve highly 
relativistic motions jjj. This model (fireball) involves 
the acceleration of matter to large Lorenz factors (7 w 
100—1000) from a compact object, and the formation of a 
highly relativistic shock. The shock converts the directed 
kinetic energy (mostly of the baryons) into internal en- 
ergy, equally distributed between electrons and protons. 
Electrons then promptly radiate their share of the in- 
ternal energy, giving rise, through synchrotron radiation 
in a suitably generated magnetic field, to the observed 
GRB. 

The relativistic environment surrounding the above- 
mentioned shock seems suitable for the acceleration of 
protons to high energies ||. In non-relativistic shocks, 
particle acceleration is a painfully slow process: parti- 
cles shuffle diffusively from downstream back to upstream 
and viceversa, each time increasing their speed infinites- 
imally (this is the modern version of the Fermi mecha- 
nism, U). Instead, in relativistic shocks, the distribu- 
tion function of non-thermal particles in the fluid frame 
is strongly collimated in the direction perpendicular to 
the shock, and they suffer deflections which differ little 
from forward/backward scattering jj]. Furthermore, at 
each cross shocking, their energies are multiplied by the 
factor 7 2 ; for the large Lorenz factors for the shocks in 
the GRBs' scenarios described above, the largest known 
in the Universe, just two or three cycles suffice to propel 
protons to energies ~ 10 20 eV. In GRBs, a large part 
of the energy loss must go through this channel, but the 
mechanism seems so fast and powerful that it has been 
suggested that the whole flux of Ultra High Energy Cos- 
mic Rays (UHECRs) observed at Earth is generated in 
GRBs |,| . 

The major rival to acceleration of UHECRs in GRBs 
is acceleration in blazars, or in the hot spots of radio 
galaxies w] . This mechanism is however hampered by the 



paucity of blazars and radio galaxies inside the Greizen- 
Zatsepin-Kuzmin limit [Q, < 100 Mpc. Thus it has 
problems explaining the rough isotropy of the directions 
of arrival of the UHECRs, and the lack of any suitable 
candidate as the site of acceleration around the direction 
of arrival of the highest energy cosmic ray observed so 
far which are not a problem for the GRB theory ||. 

However, the GRB scenario ought to find independent 
confirmation. One possible test has been proposed al- 
ready |l(J : it involves the fact that the spectrum of cos- 
mic rays emitted by a single source, as observed at the 
Earth, is wide (a power-law) if the source is a continuous 
emitter, but is much narrower if the source is explosive, 
because of energy-dependent time-delays in the propa- 
gation of cosmic rays in the magnetic field of intcrgalactic 
regions. In this Letter ; I propose a different test which 
may also reveal something about the details of the accel- 
eration process. 

Synchrotron emission by UHECRs. In models for cos- 
mological GRBs, two different scenarios for the genera- 
tion of the shock have been envisioned so far: in the first 
ipH , a shock is generated when the ejecta crash into the 
interstellar medium, much like a SuperNova. In the sec- 
ond (12), the compact object generates two expanding 
relativistic shells, endowed with slightly different Lorenz 
factors; when the second, faster shell overcomes the first 
one, collisionless shocks propagate through both shells. 
It has been argued |u| that the first mechanism may 
be responsible for bursts with smooth lightcurves, and 
the second one for spiky lightcurves. For reasons to be 
explained later, I concentrate on the first one. 

In this scenario, a total energy release Eqrb = 
-E51IO 51 erg is contaminated with a baryon mass Mb such 
that 

n ee E GRB /M b c 2 « 10 3 . (1) 
The baryons, and the Coulomb-dragged electrons, are 
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accelerated to a Lorenz factor 7 w 77, and then start a 
coasting phase which is terminated with the formation of 
a shock with the interstellar medium, of number density 
til w 1 cm -3 , for typical ISM environments. The shock 
forms approximately when a total ISM matter ps M/rj 
has been collected. This occurs at a radius 



r d = 10 18 cm El{ 3 n 



V3 -2/3 
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(2) 



At this moment, the just-formed shock splits: a forward 
shock moves into the ISM matter, with Lorenz factor 
\/2t], while a reverse shock propagates backward into the 
ejecta, causing the shell's deceleration. It can be shown 
p| that the ejecta shell thickness at this moment is given 
by r s h — Td/rj in the shell frame, where the post-shock 
magnetic field is B = 1 G n\ r}^ 1 ^ 2 , with £ parametriz- 
ing departures from exact equipartition (corresponding 
to £ = 1); for this GRB model to work properly, it is 
necessary that that £ ~ 1, which I shall assume hence- 
forth. 

It was argued in [g| that the distribution of non- 
thermal protons extends from Ei ~ rfrripC 2 10 15 eV 
to 



E u = 10 19 eV V 1/3 EU\l /6 e /2 



(3) 



all energies measured at the Earth. In the shell frame, 
the lower and upper Lorenz factors become 7/ ss 77 and 

i u « io^/io 3 )- 2 / 3 !^ 3 ^/ 6 ^ 2 . 

The spectrum of non-thermal particles behind non- 
relativistic shocks is a power-law with index p ps 2. Al- 
beit less is known about relativistic shocks, a similar con- 
clusion holds Il4|]. This is consistent with observations: 
Waxman p| has computed the injection spectrum at 
GRBs, such that, after inclusion of photopion losses j7p 
and cosmological effects, these particles fit the cosmic 
ray spectrum observed at Earth, for observed energies 
exceeding 3 x 10 18 eV. He finds good agreement for any 
index such that 1.8 < p < 2.3, consistent with produc- 
tion of UHECRs in GRBs. I shall thus take, in the fluid 
frame, dn = N i~ p dy, for 7; < 7 < j u . 

For a power-law energy distribution of non-thermal 
particles, the spectrum of synchrotron emission follows a 
power-law cx v~ a with index s = (p — l)/2, ~ 0.5 in our 
case. Thus synchrotron emission is heavily dominated 
by the high-energy end of the spectrum. Each non- 
thermal particle, in the fluid frame, emits at a typical 
frequency given by lo c — 3e 2 e£?/2m 3 c 5 , where I dropped 
the inessential dependence on the angle between the field 
and the particle mean velocity, and e is the particle en- 
ergy in the fluid frame 1=3 E/t). The high-energy end of 
the photons' spectrum is cutoff at the typical emission 
frequency lu uc of the highest energy particles, E u , which, 
scaling to the values given above, and trasforming the 
emission frequency to the Earth's reference frame, is 



E r 



1 Gevv /3 ££Nr /u £ 



2/3^-1/6^3/2 . 



for the values favoured by observations of r\ = 10 3 and 
Esi =4 m, the cutoff energy becomes E co w 300 GeV . 

While the spectral shape and the cutoff energy are easy 
to derive, the overall normalization is trickier. In non- 
relativistic shocks [jlTj the total energy emitted as radia- 
tion is of the same order of magnitude as that channeled 
into non-thermal particles; while for relativistic shocks 
no equivalent computation exists, still the arguments 
given above make it clear that relativistic shocks ought to 
be, if anything, more efficient than non-relativistic ones 
in accelerating protons. Assuming rough equipartition, 
the total flux of UHECRs at Earth is explained by the 
GRBs' scenario Thus I will equate the total energy 

released in cosmic rays Eqr to that released by the burst 
in radiation, both, of course, in the shell frame: 
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obtaining N = E G R B /rim p c 2 In 7^/7/. 

The total mass in non-thermal particles Mcr thusly 
determined is reassuringly small: we have 



P u 9 , N m v 
M C r = N m p / 7 2 df = v - 

which, using Eqs. 0, || and 7/ ps 77, gives 
M C r 1 



M b rjhiiu/ii 



< 1 



(6) 
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The synchrotron energy loss per particle per unit time 
is given by e = — 2e 4 £> 2 e 2 /(3rn 4 c 7 ). Most of the energy 
is lost by the highest energy non-thermal particles; for 
those with energy « E u , it can be shown that the ratio 
of synchrotron deceleration time to shell crossing time is 



^ = 9o.o< /3 i4 /3 r 3/2 

trr 



(8) 



which depends weakly upon all parameters except the 
efficiency £ with which equipartition magnetic fields are 
built up behind the shocks. 

The total energy loss through synchrotron emission 
by non-thermal protons is obtained by integrating the 
energy loss rate e over the particle spectral distribu- 
tion, with the normalization given above, and multiply- 
ing times the flight time across the shell thickness, which 
is the region over which the magnetic field is appreciable. 
Transforming then to the Earth frame, the total energy 
radiated is 



E sy = 3.0 x 10 49 erg 



V 



-1/3-7/6 p 5/3.3/2 
"l "51 ? 



(9) 
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For the favoured values 77 = 10 3 and E51 = 4, I find 
E sy w 3.0 x 10 49 erg. For p = 2.3, the limit of the range 
allowed by fitting the UHECRs' spectrum at Earth p5[, 
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the above computation yields E sy w 3 x 10 erg. The 
reduction of the yield in GeV photons is due to the fact 
that steep particle spectra have less energy in the high- 
energy region, where synchrotron losses are stronger. In 
summary, I expect a fraction w 0.01 of the total photon 
energy release to end up in the GeV region. 

The spectrum can now be rewritten in terms of the 
photon energy as observed at the Earth E p h in units of 
the cutoff energy, Eq. ||, x = E p h/E co , as 

dEffl = ^x~ s ok (10) 
y s 

with E sy given above. Since detectors in this energy 
range are obviously photon counting devices, it is con- 
venient to give also the particle spectrum: 

cnM = E^L x -is+D & ra 1 - 0Kl ° m x -(s+i) to 

y sE co s 

The last loose end left is to check whether the optical 
depth for pair creation against the photons emitted by 
the electrons exceeds unity: in this case, the GeV-range 
photons would be degraded in energy, and no flux near 
the photon cutoff energy might reach the Earth. I follow 
in this ref. fl§| |. The total optical depth is t 77 = r s / l /Z 77 , 
where the mean free path for photon-photon pair cre- 
ation is ^ 77 = (JTUjCph I (4m e c 2 ) 2 . Here <jt is the Thomp- 
son cross section, U 1 = Egrb / (inr^) is the photon en- 
ergy density in the shell frame, and e p h the test photon 
energy also in the shell frame. The above formula approx- 
imates the photon/photon pair creation cross-section as 
a constant, sa 3ctt/16, neglecting its decrease with in- 
creasing energies. It is thus, strictly speaking, an upper 
limit to the optical depth, which is adequate here. It can 
easily be seen with the values provided above that 

r 77 = 1.2xl0- 5 S 5 1 1 / V /3j ^V- (12) 
77 51 m e c 2 

For the usual favoured values E51 = 4,rj = 10 3 , I find 
that r 77 = 1 only for photons which, as seen from Earth, 
exceed 3 TeV. Thus there will be no energy degradation 
in situ because of photon/photon pair creation. 

This contrasts with the opposite result obtained in ref. 
Ji"8|] , but this is because they considered the other sce- 
nario for the generation of the shock, i.e., where two rel- 
ativistic shells collide with each other. In their scenario 
the shock occurs at smaller radii, so that the photon en- 
ergy density is much higher than in the scenario adopted 
here. This is the reason why I concentrated on ISM shock 
scenario. It should be stressed that UHECRs are accel- 
erated in both scenarios, the only difference being that 
no GeV photons can come out of one of them, and that, 
potentially, detection of these high energy photons may 
distinguish between the two GRB scenarios. 



Observability. I now consider the observability of the 
highest energy photons, those with energies close to E co . 
High energy photons produce pairs by collisions with 
photons of the Infra Red or Microwave Background, thus 
losing energy efficiently. This limits the range from which 
300 GeV photons can reach us to D m ss 300 Mpc M. 
The flatness of the log N — log S relation for bursts (19] 
implies that we are already seeing the edge of the GRBs' 
distribution; thus the rate of GRBs deduced from ob- 
servations, ~ 30 yr^ 1 Gpc~ 3 , is not likely to be very 
incomplete. From this, I deduce a rate of GRBs inside 
D m of 3 yr^ 1 . For these distances, the expected fluence 
from Eq. [| is 3 x 10~ 6 erg cm~ 2 . 

There are currently no experiments which can detect 
showers initiated by primaries with energies around E co 
with the required sensitivity, but the next generation of 
high-altitude (> 4000 m a.s.l.) detectors currently be- 
ing planned will have relevant detection thresholds of 
sa 10 GeV. Since I argued that the spectrum is ex- 
pected to be very flat, most energy, though not most 
counts, will be deposited above this detection threshold, 
so that the limiting factor will be the experiments' detec- 
tion surface. As an example, ARGO can detect sig- 
nals down to flux levels of *=s 3 x 10 -7 erg cm~ 2 for spec- 
tra extending out to 300 GeV, provided the bursts last 
1 s, and the spectra are flatter than s sa 1.5. However, 
with this low detection threshold, a burst with fluence 
3 x 10~ 6 erg cmT 2 could be detected over background 
noise even if it lasted 100 s. 

Implications. High energy photons will be emitted as 
long as protons are accelerated to very high energies, 
which only requires the presence of large magnetic fields. 
There is no obvious reason why the field should decay 
on the timescale of the burst. In fact, first, as seen from 
Earth, the shell remains relativistic for about a month af- 
ter the burst uM , which implies that relativistic electrons 
will be available in the post-shock region to generate a 
non-negligible magnetic field. Second, after the burst, 
the electrons exchange energy with protons through a 
variety of processes. So long as electrons are thusly kept 
relatively hot, their chaotic motion may maintain an ap- 
preciable magnetic field. Since processes coupling elec- 
trons to protons are relatively inefficient in transfering 
energy, some magnetic field may be left for some time 
after the burst. So observations of the secondary burst 
discussed in this paper may reveal burst durations ap- 
preciably longer than those of their lower energy coun- 
terparts. 

Another important issue is whether this emission may 
be masked by Inverse Compton effects: it was shown E| 
that photons with energies as high as 10 (n/100) 6 GeV, 
very incertain because of the steep dependence upon rj, 
can be produced this way. In this case, the expected 
spectrum is simply that of lower energy, seed phtons, 
which get a kick to higher energies: typically the photon 
number spectrum is cx v~ 2 , different from Eq. nTX This 
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furthermore implies that the energy released per photon 
energy decade is constant, and thus exceeds the estimate 
of Eq. ^| by about two orders of magnitude. Thus, spec- 
tral steepness and fluence allow an easy discrimination 
between the two mechanisms. 

The above argument finds an interesting application at 
intermediate energies, rs 1 GeV. One may in fact won- 
der whether the low energy tail of the emission derived in 
Eq. |l^ has been observed already, since some bursts do 
show detectable fluxes at these energies . This can be 
certainly excluded for GB 930131, which has comparable 
fluences in the low-energy BATSE spectrum (w 1 MeV) 
and in the GeV region, with nearly identical spectral 
slopes 23 1 . The situation is however less clearcut for 



other bursts, like GB 940217, GB 910503, which have 
fluences, at their respective highest energies, much lower 
than those emitted at lower energies, and thus in keeping 
with the estimate of Eq. [|. It should also be noticed that 
all these bursts show, at the highest energies, burst du- 
ration much longer than at low energies, a fact naturally 
accounted for in this model, as discussed above. 

A mechanism for the production of delayed GeV/TeV 
photons from GRBs based upon the presence of UHECRs 
has been proposed [Q, but it differs greatly from the 
present one: there UHECRs produce photons through 
photopion processes off CMB photons in flight from the 
site of the burst to our detectors, while in this model the 
emission mechanism is synchrotron in situ. 

Summary. Physical conditions at the relativistic 
shocks which give origin to GRBs are surely favorable 
to acceleration of high energy particles. If the energy 
in UHECRs is comparable to that released in photons 
(and thus if UHECRs from GRBs account for the whole 
flux of UHECRs at Earth) , the fluxes are those predicted 
above. Otherwise, the observed flux of GeV photons will 
allow measurement of the fraction of energy channeled 
into non-thermal particles. 

I have shown that the synchrotron emission spectrum 
from protons extends to E co ps 300 GeV (Eq. [|) and 
should be rather flat (Eq. [Io|), that the total energy re- 
lease is E sy ss O.OIEqrb (Eq. ||), and that a few such 
events per year will surely become observable with the 
next generation of high altitude, small size air shower 
detectors. I have also argued that all GeV emission ob- 
served to date from GRBs is compatible with this model. 
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